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We propose the another, in principe nonperturbative, method of the evaluation of the Wiener 
functional integral for ip'^ term in the action. All infinite summations in the results are proven to be 
convergent. We find the "generalized" Gelfand-Yaglom differential equation implying the functional 
integral in the continuum limit. 



We follow the definition of a functional integral as a limit of a sequence of the finite dimensional integrals. This 
method has no problems with the definition of integral measure, because for the finite dimensional integrals the 
integral measure is defined correctly. Calculating the finite dimensional integrals we must solve the problem of the 
calculation of one dimensional integrals 

/r=/ ^.exp{-(a.^ + 5.^ + c.)}, (1) 

— oo 

where Rea > 0. The conventional perturbative approach rely on a Taylor's decomposition of x'^ term with consecutive 
replacements of the integration and summation order. These integrals can be calculated, but the sum is divergent. 
However, Ji = /i(a, 5, c) is an entire function for any complex values of b and c, since there exist all integrals 



dl'dj^'hia, b, c) = (-1)"+" J dx exp{-(aa;'* + bx' + cx)} . 

— oo 

Consequently, the power expansions of Ii = /i(a, 6, c) in c and/or 6 has an infinite radius of convergence (and in 
particular they are uniformly convergent on any compact set of values of c and/or b). In what follows we shall use 
the power expansion in c: 

Ii = Y ^-i- / dx x" exp{-(ax4 + bx^)} . (2) 

The integrals appearing here can be expressed in terms of the parabolic cylinder function D^{z), v — —m— 1/2, (see, 
for instance, For n odd, due to symmetry of the integrand the integrals are zero while for even n = 2m we use: 

D-,n-i/2{z) = ^^^^1/2) i, '^'^ ' exp{--a;2 - zx} . (3) 

Explicitly, for the Eq.(|21) we have: 

^ m— ^ 
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we have used the abbreviation: 



The sum Q is convergent for any values of c, h and a positive 

Using such expansions we shaU calculate the continuum unconditional measure Wiener functional integral: 

[X>(p(x)]exp(-5) , 

corresponding to the continuum action with the fourth order term: 



c/2 



I dr 



(5) 



Following the standard procedure, we divide the integration interval into N equal slices. We define the N— dimensional 
integral by the relation 0: 




N 



exp<^-^A 



c/2 



A 



(6) 



where A = P/N . The continuum Wiener unconditional measure functional integral is defined by the limit: 



Z = lim Zjv 



Applying to © the relation (|3Jl, we obtain, 0: 



Zat = [27r(l + 6AVc)] ^ [2^(1/2 + 6AVc)] ^'^Sn 



with 



N 



E n 

k\,--- .kN -1—0 i— 



{2ki) 



■r(fc,_i + h + i/2)p_fc,_^_fc^_i/2 (z) 



(7) 



(8) 



where ko = kN = 0, p = (l + 6A^/c) ^, z = c{l + bA^/c)/^/2aA^ and the argument of the last term I'-fc„_i_i/2(-Z;ast) 
is ziast = c(l/2 + 6A^/c)/ ^/2aA^. It is worth while to stress that the factor p is independent of the coupling constant. 
The coupling constant dependence enters into this formula only through the argument z of I? functions. 

The question of the convergence of the Eq.® is important. Let a^^ is the ki dependent part of the argument of 
the product in Eq.JSJl: 



ip) 



2ki 



{2ki) 



■r(fc,_i + h + i/2)i?_fe,_,_fc,_i/2 (z)r(fc,+i + fc, + ii2)v_k,+^-k,-i/2 (z) 



(9) 



Using the asymptotic for parabolic cylinder functions and Stirling formula we prove that au^ approaches to zero in 
the asymptotic region of ki as 

^Dpk, 



ki\ exp (vfci) 

where D and F are finite numbers P|. This asymptotic of a^. is sufficient for a proof of the uniform convergence 
of the series not only for single ki in JH)), but for arbitrary tuple {ki} of indices. By the same method we can prove 
the uniform convergence of the summation over two, three, iV — 1 summation indices ki in the equation ©. This 
important conclusion indicates, that in the formula for Zj^ the [N — l)-tuple summation over kiS is convergent. 
Following the idea of Gelfand and Yaglom the functional integral in the continuum limit is defined by the relation 

hm Zm = — , , 
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where F{(3), for our case, is the solution of the differential equation 

3"^ 8 3 /2h 3'^ \ 

—Fir)+A-Fir)-lnS{r)^Fir)[^-^2^lnSir)J , (10) 

calculated at the point /? (the upper limit of the time interval in the action (jSJ)- Ea. pU|l has to be supplemented by 
the initial conditions: F{0) — 1 and 

dF{T)/3T\r=0 = . 

The function S{t) is given as the limit N ^ oo oi Sn given by 0. 

The result of the summation JHJl (replacing the N dimensional integration) is an exact relation, calculated without 
any approximation. The multiple summations suppress this advantage somewhat, therefore we shall discuss the ki 
summations in the formula • The details of this calculations are presented in . 

These summations can be provided by the help of the asymptotic expansions of the parabolic cylinder functions 0] : 

J (k- + 1/2) 

{J is fixed finite number, z is large), and the summation relation for the generating function 01 : 



( \ 

e^^" E ^ D^^^uix) = D^^ {x - t) . (12) 

fe=0 

The symbols {ki + l/2)„. and iv), are the Pochhamer's numbers defined by: 

{u)^^T{u + k)/nk) . 

In the summation over single fcj in © , corresponding to the summation over . in (jS)) , there appear a product of two 
D functions possessing as the function's subscript the summation index. We propose to use the asymptotic expansion 
for one of them. The asymptotic formula can be used in summation over index ki provided that the summation runs 
over finite number of the terms, and the index of the parabolic cylinder function, and its argument obey the relation: 

k, < iVo <| z I , (13) 

A^o is an upper bound of summation variables in ((SJ. 

This requirement for an application the asymptotic formula is violated if the summation index approaches infinity 
in We can overcome this problem thank to the uniform convergence of the summations in ^ proven in [j]. 
We truncate the summation over ki in © but we demand that the partial sum is close to the full infinite one. To 
minimize the effect of the truncation on the final result in the continuum limit, we have to choose the maximal value 
of the index ki such that the remainder of the sum (jSJ satisfies the inequality 

oc 

I E I < ^(^) ' (14) 

k,=No + l 

with e{N) approaching zero for large N so that N^e{N) in the limit — > oo. We have chosen the third power 
of N in the inequality in order to eliminate the influence of remainders (|14|) in Gelfand-Yaglom procedure for the 
calculation of continuum limit. We estimated the remainder (|14|) and shown that both conditions (|13|l and H14|l for the 
application of the proposed summation method are valid simultaneously. In this estimate we utilized the asymptotic 
form of the parabolic cylinder functions with double asymptotic properties proposed by N. Temme 5]. 



The result of summations gives for © the following relation 



N 



-1/2 



2n = {11 [2(l + 6AVc)a;.] 



J 

E 



(-1)^ 



(15) 
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with the symbols {N)2j j satisfying the following recurrence relation: 



M / \ 2A 



A— i—max[0, p — 2(/^— A)] 



a: 



^'^ . _ ^ A (1/2), 







"-l)2,.p = «^ = ( ^ )77^, (16) 



(1/2) 



where 



^2 

w,: = 1 



2 2(l + 6A2/c)' 



^ i-2(TtW) • ^'^^ 

-Zat in p5|) is the N ~ th approximation of the functional integral. The continuum limit of Zn we calculate by the 
procedure proposed by Gelfand- Yaglom following the analogous steps as for the calculation of functional integral for 
the harmonic oscillator 0|. This means, that for finite N we construct the difference equation for p5|l . which in the 
continuum limit A'^ ^ oo is reduced to differential equation. The method and the solutions are discussed in 0, here 
we present the result only. 

The function S{t) entering differential equation (|10|l is calculated by: 

Sir)-l^^XlAW^^^''''^''^'^^^^ (18) 

with the expansion factor 

1 (0'P\^ 

(2i2A3)^ ^ V"?J 

finite in such limit. The additional power A"^^ certifies that in the relation for (A^)2/i, o only the leading terms in 
powers of 1/A survive the continuum limit. 

Equation (|10|l can be simplified by the substitution 



F{r) ^(") 



For t/(r) we find the equation 



— y(r) = y(r) 



(19) 



accompanied by the initial conditions 



y(t>) = s"(0). (20) 

The evaluation of function S{t) is not yet finished. We present the result, related to the usual perturbative 
expansion, to the lowest order terms in power expansion in the variable a. Up to the term linear in coupling constant 
a we have found 0: 

5W(r)= hm^^ (A^)2,,o = l-^^ {tanh(r7) + r7[3tanh2(r7)-l]} , 
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where 7 = yj2b/c. 

We see, that for 6 > in the harmonic oscillator limit, (i.e. a = implying S'^^^t) = 1), Eqs. Ipll CHI reduce to 
Gelfand- Yaglom equation for the harmonic oscillator. We find that in lowest nontrivial order in a the Eq. p9|l is of 
the modified Bessel type one. Then the functional integral can be expressed in term of the modified Bessel functions. 

Inserting the asymptotical perturbative expansion S*^^^ (r) into equation (|19() . we calculate functional integral beyond 
simple perturbative expansion. We can render our result as a partial resumation of the simple perturbative series. 
Moreover, this procedure allows the parameter h take both positive (anharmonic oscillator case) and negative (Higgs 
case). A similar procedure is familiar in the renormalization group calculations proposed by Gell-Mann and Low. In 
Callan - Symanzik equation for running coupling constant the [3 function is calculated perturbatively, however, the 
equation is solved non-perturbatively. 
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